Minimum free-energy path of homogenous nucleation from the phase-field equation 
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The minimum free-energy path (MFEP) is the most probable route of the nucleation process 
on the multidimensional free-energy surface. In this study, the phase-field equation is used as a 
mathematical tool to deduce the minimum free-energy path (MFEP) of homogeneous nucleation. 
We use a simple square-gradient free-energy functional with a quartic local free-energy function as 
an example and study the time evolution of a single nucleus placed within a metastable environment. 
The time integration of the phase-field equation is performed using the numerically efficient cell- 
dynamics method. By monitoring the evolution of the size of the nucleus and the free energy of the 
system simultaneously, we can easily deduce the free-energy barrier as a function of the size of the 
sub- and the super-critical nucleus along the MFEP. 

PACS numbers: 64.60.Q-,64.60.qe 
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I. INTRODUCTION 

Nucleation is a very basic and ubiquitous phenomena 
that occurs in phase transformation of various materials^. 
A small embryo of new stable material appears from the 
fluctuation within a metastable material, and it grows 
into a small cluster called a nucleus. Then it overcomes 
the free-energy barrier and will grow indefinitely. The 
basic process of the formation of a single nucleus called 
nucleation is difficult to study even theoretically as the 
process involves non-equilibrium and transient states. In 
much of the literature^i^i^ only the critical nucleus that 
corresponds to the cluster at the top of the energy barrier 
has been considered and used to deduce the activation 
energy of the nucleation rate. 

The whole landscape of the free-energy surface for a 
single nucleus is impossible to study because of the huge 
degree of freedom even for a small cluster. Intuitively, 
however, it is believed that the sub-critical nucleus as- 
cends the free-energy surface along the valley to the top 
of the barrier. The top of the barrier is in fact the saddle 
point of the free-energy surface. The nucleus at the sad- 
dle point is called the critical nucleus. It surpasses the 
barrier at the saddle point and descends along the valley 
as a super-critical nucleus and will grow indefinitely. This 
minimum free-energy path (MFEP) is the most probable 
route of the reaction for the nucleation even though the 
real path might wander around this MFEP through a 
thermal fiuctuation or some other effect. 

In order to extract information about the free-energy 
surface for the sub- and the super-critical nucleus, Weak- 
liem and Reiss^ proposed a theoretical model called the 
molecular theory of nucleation. In this theory, the {N, A) 
cluster model where a fixed number of molecules N con- 
fined within a container of the radius A (or the volume 
V) are considered. Using a physically plausible boundary 
condition one can study the projection of the free-energy 
surface onto the two-dimensional {N, A) (or {N, V)) space 
which is expected to capture the essential features of 
the real free-energy surface in multi-dimensional space. 



This theory was originally formulated using the Monte 
Carlo method and was reformulated using the density- 
functional theory (DFT)^i. However, it is not obvious 
how to restore this projected free-energy surface in the 
two-dimensional (iV, A) space to the true free-energy sur- 
face and barrier in the multi-dimensional space. There 
are also many problems inherent to the confined system^. 
Another theoretical model to use the pressure P instead 
of the virtual radius A has also been studied^. 



The simplest way to study this free-energy surface is to 
look at only the MFEP for a single nucleus as the function 
of the number of molecules N asa. reaction coordinate. In 
this case, the fluctuation of the size A or the volume V for 
the given number of molecules N considered in the {N, A) 
cluster model of the molecular theory^ii is neglected. 
The Monte Carlo simulation with biased samplingiSiiii 
can be used to trace the MFEP. Similarly, by choosing 
the parameterized spherical density proflle, or by using 
constraints in the minimization process we can calculate 
the MFEP as the function of the number of molecules 
using DFT. 



In this paper, we will propose an alternative novel 
method to deduce the MFEP using a much simpler phase- 
field equationi^ for the non-conserved order parameter. 
Since the time integration of the highly non-linear phase- 
field equation requires a certain amount of computational 
resourcesi^, we will further simplify the equation using 
the cell-dynamics methodic. The format of this paper is 
as follows: Section |TT] is a brief review of the phase-field 
equation and the free-energy model used to study the 
MFEP. In section [1111 the cell-dynamics method is used 
to integrate the phase-field equation and the MFEP is 
extracted from the time evolution of the size and the free 
energy of a single nucleus. We will conclude this paper 
with our final comments in section [TVl 
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II. PHASE-FIELD EQUATION AND THE 
MINIMUM FREE-ENERGY PATH 



In order to study the evolution or the regression of 
a single nucleu a^^i^^ and the reaction-path of nuclcation 
called the minimum free-energy path (MFEP). wc will 
use the partial differential equation called the phase-field 
equation for the non-conserved order parameter: 



'~dt 



SF_ 



(1) 



where ip is the order parameter called phase field. F is 
the free-energy functional (grand potential) usually ex- 
pressed by the square-gradient fornt^ '^^i^^ , and r speci- 
fies the time scale of evolution. This phase-field equa- 
tion is attractive as the kinetics of phase transformation 
is driven by the relative stability of each phase and by 
the topology of the associated free-energy surface F. In 
particular, the special solution of Eq. ([T]) is the interface- 
controlled growth of a spherical or a circular grain with 
nearly constant velocityi^ in d = 3 and d = 2 dimensions 
respectively. The dynamics driven by Eq. ([T]) is also at- 
tractive to study MFEP as it always guarantees that the 
total free energy decreases monotonicallyi^: 



dF 
'dt 



SF dip 



(2) 



Mathematically, the evolution of the nucleation is the tra- 
jectory of the dynamical system in infinite dimensions de- 
scribed by Eq. ([T]) in the phase space whose Lyapunov 
functional is the free energy F. 

Since Eq. ^ implies the steepest descent of the trajec- 
tory of "0 which represents a single spherical or a circular 
nucleus on the free-energy surface F, the phase-field ki- 
netics equation Eq. ([T]) for the non-conserved system can 
be used as a mathematical tool to study the MFEP of nu- 
cleation not only for the non- conserved order parameter— 
but also for the conserved order parameter—. Equation 
([T]) can also be interpreted as an over-damped dynamics 
on the potential surface F in the functional space ip. In 
fact, our use of Eq. ^ to explore the MFEP is analogous 
to the over damped Langevin dynamics used in the string 
method proposed by Qiu et al.— to study the transition 
pathway of nucleation in the capillary condensation. 

For the conserved order parameter, however, the ki- 
netic equation 



(SF 



(3) 



should be used instead of Eq. ([T]). Also, in the more rig- 
orously formulated dynamical density-functional theory 
(dynamical DFT) ^^i^^ , the kinetic equation is given by 



dp _ / „5F 

it - ^ r ^ 



(4) 



where the order parameter p represents the density. 
These kinetic equations Eqs. ^ and Q could also be 



used to search for the free-energy surface as they can also 
satisfy the inequality similar to Eq. Q which guarantees 
a monotonically decreasing free energy as a function of 
time. 

To be exact, the real physics of the nuclcation in a con- 
served system will be described by Eqs. ([3]) or ([5]). How- 
ever, these equations have a shortcoming in that the dy- 
namics will be complex due to the depletion effect as the 
integrated order parameter or the total density must be 
conserved. Also, the mathematical structure of Eqs. ^ 
or (HI implies that they do not trace the free-energy sur- 
face F itself in a straight forward maner. In contrast, the 
phase-field equation ^ can trace the free-energy surface 
and the MFEP directly. Therefore, we will use Eq. ([T|) 
to study the MFEP even if it may not have a physical 
meaning for the nuclcation in conserved systems such as 
liquid-vapor nuclcation. 

As an example, wc will employ the square-gradient 
model for the free-energy functional: 



FM = 



/(V') 



dr 



(5) 



where the local part of the free energy /(ip) proposed by 
Jou and Lusk^^ is used: 

/(^) = i^^^(0-l)^ + |(^-^), (6) 

which represents the two-phase system with two phases 
characterized by ip = and ip = ^ shown in Fig. [TJ The 
parameter e controls the relative stability of one phase at 
■0 = 1 with the grand potential f{ip = 1) = — e/4 relative 
to another phase at -0 = with the grand potential f{ip = 
0) = 0. Wc will only consider the case when e > 0. Then 
the phase at ipm = is the metastable and the one at 
-03 = 1 is stable. The free-energy difference A/ (Fig. [1]) 
is given by 



A/ = /(0)-/(l) 



(7) 



(8) 



The position ipi, of the free-energy barrier is given by 

1 3e 

and its height (Fig. [T|) is given by 

(77 - 3e)' (e + 77) 



fb = f ii'b) 



(9) 



with /ft 77/64 when e <C 77. Then the spinodal point 
defined by /& = is at 



e = 77/3, 



(10) 



and the binodal is at e = 0. A typical shape of the free- 
energy function f{tp) near the binodal is shown in Fig.[T] 
Now, the special solution of the phase-field equation 
Eq. ([T]) for the traveling wave of the form -0 (A) with 
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FIG. 1: The local part of the free energy Eq. when 77 = 0.7 
and e — 1/75 ~ 0.01333. We consider e > so that the phase 
at i/) = 1 is always stable. 



Equation (|lip can also be integrated directly 



dX [2dX^ 



dij 
dX 



to give another expression for the velocity v 

DAf 

V = 

TO" 

expressed by the surface tension a defined by 



dX 
(18) 

(19) 




^4-\dx. (20) 
dx J 



The direct integration of Eq. (|20p with the free energy 
given by Eq. ([6|) is analytically possible only at the two- 
phase coexistence e = 0, and the surface tension a is 
gives 



X = r — R{t) having a single spherical or a circular shape 
of the radius R{t) satisfies the differential equatio n^^i^^ 



D 



d^ijj 
dJO 



where 



) — 
dX 



dR 



f ii') - 



D- 



R 



(11) 



(12) 



and A = 0,1,2 for d = 1,2,3 dimensions. The nucleus 
grows or shrinks according to 



dR 



1 



Rc 
R 



(13) 



from Eq. Now the nucleus will grow if the radius 

R becomes larger than the dynamical critical radius Rc 
given by 



Rr = 



DA 



(14) 



otherwise it will shrinki^. The time evolution of the ra- 
dius R{t) is given hy^ 



vt = (i?(t) - i?o) + Rc In 



R{t) - Rc 
Rq — Rc 



(15) 



from Eq. p^ . where Rq is the initial radius at i = 0. 
For a sufficiently large nucleus {R{t) 00), the radius R 
grows with a constant velocity dR/dt = v given bjJ^iS^ii^ 

(16) 

and the order parameter profile ipi^) the form 

1 




^(X) 



exp 



'^X 

D ^ 



(17) 



(21) 

Inserting Eqs. ([7]) and pT|) into Eq. p9|) . we can recover 
Eq. (fT6|l . Now we will simulate the evolution of a sin- 
gle nucleus using the phase-field equation Eq. (jT]). By 
monitoring the decreasing free energy of a single nucleus 
as a function of its radius, we can trace the MFEP for 
the shrinking sub-critical and the growing super-critical 
nucleus. 

In the classical nucleation theory (CNT) , the minimum 
free-energy path (MFEP) is given by the work of for- 
mation WcNT of the nucleus calculated by assuming a 
spherical or a circular shape and a sharp interface with 
the size-independent surface tension a: 

VFcNT = -yi?'A/ + 4^i?2^, {d^3), 

= -nR^Af + 27rRa, (d = 2), (22) 

for (spherical) d = 3 and (circular) d = 2 dimensions. 
This is the "minimum" free-energy path within the CNT 
as we have assumed a spherical or a circular shape of 
the minimum surface area. Other shapes with a larger 
surface area and the same volume naturally give a higher 
free energy as expected from the second term of Eq. (|22p . 
By maximizing this free energy Eq. (|22p we can deduce 
the thermodynamic critical radius Rc that corresponds 
to the saddle point of the free-energy surface for a single 
nucleus: 



(23) 



where the meaning of A is the same as that in Eq. (|12|) . 
Then the free-energy barrier (activation energy) of nu- 
cleation rate is the maximum of the MFEP given by 
Eq. dllD: 



W, 



16tt 



CNT 



3 (Af) 
\ id ^2) 



(d = 3). 



A/ 



(24) 
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The thermodynamic critical radius of the CNT in 
Eq. (|23p is the same as the dynamical critical radius in 
Eq. (fT4|) since both expressions reduce to 



Rc = 



V2^ 
3e 



(25) 



for A = 1 (d = 2 dimension), for example, from Eqs. jT]), 

mi and im. 



III. 



RESULTS OF NUMERICAL SIMULATION 
USING CELL DYNAMICS 

A. Cell-dynamics simulation 



Since Eq. ([T]) is a highly non-linear partial differential 
equation, we will not attempt to solve it directly. Instead, 
we follow Oono and Pur i^^i^^ and transform this partial 
differential equation into the space-time discretized dif- 
ference equation called the cell-dynamics equation. Their 
transformation docs not mean an accurate numerical ap- 
proximation to the original phase- field equation Eq. ([T]). 
Rather, they originally aimed at simulating the global 
picture of the kinetics of the phase transformation gov- 
erned by the phase-field equation within the framework 
of the discrete cellular automata. However, later workers 
have found that this difference equation can successfully 
reproduce the time evolution of the original phase-field 
equation driven by a subtle balance of different phases 
represented by the free-energy surface Jis, 24^27 ^ 

According to the cell-dynamics method, the partial dif- 
ferential equation Eq. ([1]) is transformed into the finite 
difference equation in space and time of the following 
form 



ij{t + l,n) = M[i;{t,n)] 



(26) 



where the time t is discrete integer and the space is also 
discrete and is expressed by the site index (integer) n. 
We have eliminated the time scale r because the cell- 
dynamics equation Eq. (j26p is a coarse-grained approx- 
imation to Eq. ID) in time and space so that the time 
scale is irrelevant. The mapping M is given by 

M[ij{t, n)] = g{iP{t, n)) +D[« tP{t, n) » -^{t, n)] 

(27) 

where the definition of << * >> for the two-dimensional 
square grid is given by 

« V^(t,n) »= i ^ V^(i,i) + i^ ^{t.i) (28) 
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with "nn" means the nearest neighbors and "nnn" the 
next-nearest neighbors of the square grid. We will only 
consider the nucleation in d = 2 dimension. 

Instead of the original hyper-tangential map function 
5 (V') = ip — 1.3 tanh-i/ i^^i^^ , we use the map function g 
that is directly related to the free energy f{ip) through 



(29) 



This replacement is essential to study the nucleation and 
the growth when a subtle balance of the relative stability 
of two phases realized by the topology of the free-energy 
surface / plays a crucial role. 

Ren et ali^ argued that by using the map function 
like Eq. ((29|) instead of the original hyper-tangential 
map function one can easily include the effect of asym- 
metry in the free energy function / and, hence, can 
take into account the asymmetric character of the two 
phases within the framework of the cell dynamics. Sub- 
sequently, the present author has demonstrated that by 
using the map function Eq. derived from the free- 
energy function / (ip) one can simulate the growth of 
a single circular nucleus to confirm the analytical for- 
mula— for the growth velocity^l. Furthermore, by using 
the map function in Eq. (j29p . the nucleation and the 
growth process of multiple nuclei can also be success- 
fully simulated not only in the site saturation rcgimo^l 
but also in the continuous nucleation regime>i^ to confirm 
the KJMA (Kolmogorov-Johnson-Mehl-Avrami) picture 
of phase transformation^^. It is now well established that 
this cell-dynamics method can serve as a simple integra- 
tion scheme to study the kinetics of phase-field model 
even though the method is not guaranteed^^ to be an ac- 
curate numerical approximation to the original equation 
Eq. HD. 

Since the direct numerical integration of the partial 
differential equation Eq. ([1]) requires a finely tuned inte- 
gration schem e^^'^^ and a significant amount of compu- 
tational resourcesi^, we will use this cell-dynamics equa- 
tion Eq. (|26p instead of the original phase-field equation 
Eq. ID) to study the evolution of a single nucleus in this 
paper. 



B. Numerical results 

Initially a circular nucleus with the order parameter 
-0 = 1 and the radius Ri is placed within the metastable 
environment with V' = in d = 2 dimensional space. 
Then the cell-dynamics equation Ec^. (P5|) is solved to 
simulate the evolution of the initial nucleus. If the initial 
radius Ri is larger than the dynamical critical radius Rc, 
the nucleus is expected to grow indefinitely otherwise it is 
expected to shrink and disappear. This critical radius Rc 
is also the thermodynamic critical radius of the critical 
nucleus. Wc will use three sets of the free-energy param- 
eters shown in Table HI The parameter D in Eq. ^ is set 
to 1/2 throughout this paper. We use small e to make 
the critical radius Eq. |25|) large. Therefore we will only 
consider the thermodynamics state closer to the binodal 
than the spinodal where e = 77/3 must be large. 

Figure [2] shows the cross sections of the growing super- 
critical nucleus when 1] = 1.0 and e = 0.01333. Initially 
a circular nucleus of ^jj=l with a radius i? = 18 and a 
sharp interface is prepared. Its radius increases and the 
interface become diffuse in the course of evolution. The 
nucleus with initial radius R = 14—17 neither grows nor 
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FIG. 2: Cross sections of the growing super- critical nucleus 
when rj — 1.0 and e — 0.01333. Initially a circular nucleus of 
^p=l with a radius R — 18 and a sharp interface is prepared 
at the center of 256x256 cells. 



TABLE I: Three sets of the free-energy parameters rj and e 
used together with the critical radius Rc and the interfacial 
velocity v from theoretical predictions by Eqs. (|25p and (|16p 
and from our cell-dynamics simulations. 



V 


1.0 


0.7 


1.0 


e 


0.01333 


0.01333 


0.02222 


Rc 


25 


21 


15 


Rc (Simulation) 


14-17 


13-14 


8-10 


V 


0.020/r 


0.024/r 


0.033/r 


V (Simulation) 


0.0070 


0.0098 


0.0121 



shrinks, while the nucleus with initial radius smaller than 
i? = 13 always shrinks. The critical radii i?c = 14 — 17 in 
Table [I] estimated from our cell-dynamics simulation are 
those radii for which the initial nucleus neither grows nor 
shrinks. The critical radius are slightly smaller than the 
theoretical prediction from Eq. tabulated in table I 
that is probably due to the discrete difference Eq. pS)) 
used in the cell-dynamics method. 

The time evolution of the radius R of the nucleus de- 
fined as the distance from the center of the nucleus to 
the position of the cell with i^iR) — 0.5 shows an almost 
straight line as a function of time expected from the so- 
lution of the differential equation Eq. ([15]) for the super- 
critical growing nucleus {R > Rc) as shown in Fig. [31Ja). 
The velocity v deduced from the cell-dynamics simula- 
tion cannot be compared with the theoretical prediction 
from Eq. as the time scale t in Eq. ([T]) is unknown. 
The same problem of the time scale has already been no- 
ticed by the author—. However, we observe from Table |T] 
that the time scale t is almost constant and is around 
T = 2.5 -2.8. 

In contrast to the super-critical nucleus, the radius of 
the shrinking sub-critical nucleus shown in Fig. [3]Jb) does 
not decrease linearly in time. The radius remains con- 
stant for some periods, then it starts to decrease grad- 
ually. As the radius is getting smaller and smaller it 
decreases more rapidly. This non-linear behavior is due 
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(a)Time evolution of the radius of the 
super-critical nucleus. 
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(b)Time evolution of the sub-critical nucleus. 



FIG. 3: (a) The time evolution of the radius of the grow- 
ing super-critical nucleus for the three sets of potential pa- 
rameters in Table U] (b) Time evolution of the radius of the 
shrinking sub-critical nucleus. 



to the capillary pressure: The nucleus will be pushed in- 
ward by the curvature correction to the pressure. Since 
we use the discrete cell dynamics, we also observe the 
staircase structure in Fig.[3]Jb). 

Similarly, the time evolution of the total free energy 
given by Eq. ([5]) of the sub- and the super-critical nu- 
cleus can be calculated using the temporal order param- 
eter profile ipit, n) = ip{x,y) where n = (x,y) is the site 
index in two dimension. The necessary gradients are ap- 
proximated by the finite differences 



dx 
dy 



ipjx + l,y)- ipjx - l,y) 
2 

tlj{x,y + 1) - tlj{x,y - 1) 



(30) 



Figure [4] shows the time evolution of the total free en- 
ergy F defined by Eq. (O as the function of the time t. 
The free energy of the growing super-critical nucleus in 
Fig. mja) and that of the shrinking sub-critical nucleus 
in Fig. \Mi>) decrease monotonically as the functions of 
time according to the prediction of Eq. ^ except at the 
early stage of growth in Fig. [UJa). A small increase in 
the free energy of the super-critical nucleus is probably 
due to the numerical errors of the coarse-grained deriva- 
tive in Eq. (PD)) . The initial nucleus at i = has a much 
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high energy which is indicated by an isolated symbol at 
i = in Fig. [5] because it has an artificial shape with a 
step- function interface which gives an unrealistically high 
surface tension from Eq. (|20p . 
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(a)Time evolution of the free energy for the 
super-critical nucleus. 



S 69 



Sub-critical nucleus 



1^=0.10 8=0.0133 

oooooooooooooooooooot^ 
ri=0.07 6=0.0133 

ti=0.10E=0.0222 




i^ooooooc 



500 1000 1500 2000 2500 3000 3500 4000 
Time t 

(b)Time evolution of the free energy for the 
sub-critical nucleus. 



FIG. 4: (a) The free energy F of the growing super-critical 
nucleus as the function of time t. (b) The free energy of the 
shrinking sub-critical nucleus. Both these free energies de- 
crease monotonically in accordance with the theoretical pre- 
diction in Eq. Q except at an early stage of growth. 

Combining Figs. [3] and [H we can deduce the total free 
energy F of nucleus as the function of its radius R shown 
in Fig. m Since the locus of the phase-field equation 
Eq. ([T|) in the phase space always follows the route to 
lower the free energy of a single nucleus, the curve in 
Fig. [5] represents the minimum free- energy path (MFEP) 
of nucleation. We have also plotted the MFEP calculated 
from Eq. ((^H) of the classical nucleation theory (CNT) in 
Fig. O The overall shapes of the MFEP from our simula- 
tions and those from the CNT are similar. In contrast to 
the dynamical critical radius Rc shown in Table |T1 the po- 
sition of the thermodynamic critical radius Rc estimated 
from Fig. [5] of our simulation is almost the same as that 
from the theoretical predictions Eq. ((25)) based on the 
CNT. 

However, the absolute magnitude of the free energy 
from our simulation is always higher than that from the 
CNT prediction. It is well recognized that the CNT 
prediction is incorrect near the spinodal where the free- 



energy barrier should vanish while the barrier from the 
CNT erroneously remains finite^'^'^. Since the phase-field 
model is correct near the spinodal, the free energy from 
our simulation is expected to be lower than the CNT pre- 
diction at least near the spinodal. Since we are closer to 
the binodal than the spinodal, the free-energy barriers 
from DFT could be higher than the CNT predictions. In 
fact, many other theoretical approaches seem to predict 
the same trend as ours near the binodal. For example, 
Oxtoby and Evans^ predicted that the free energy calcu- 
lated from their non-classical DFT for the Yukawa fiuid 
can be higher than the CNT prediction near the bin- 
odal. Within the framework of the square-gradient model 
for the liquid-vapor nucleation^, the present author has 
demonstrated analytically using a double-parabolic free 
energy that the asymmetry of the free energy between the 
gas and the liquid phases accounts for the asymmetric be- 
havior of the effective surface tension of the nucleus and, 
therefore, of the free energy barrier. In this liquid- vapor 
nucleation model^, the free energy barrier for the vapor 
bubble is always lower than the CNT prediction while 
the barrier for the liquid droplet can be higher than the 
CNT prediction near the binodal. We anticipate that a 
similar asymmetry in our model free energy Eq. ([6|) will 
account for the free energy barrier higher than the CNT 
prediction in Fig. [51 



10 



A Tl=0.10 

£=0.0133 




20 30 40 

Radius R 



FIG. 5: The free-energy barrier along the minimum free- 
energy path (MFEP) of nucleation deduced from Figs. [3] and 
m The solid lines are the CNT predictions from Eq. (ggj). 



Despite the ambiguity in the definition of the radius 
R of the diffuse interface in the phase-field model, the 
overall shape of the free-energy barrier in Fig. [5] is sim- 
ilar to the one from the CNT as well as to the results 
of other researchersiSiiiii^. Naturally, no anomalous in- 
stability or catastrophe which is observed in the (N, A) 
cluster model^ for the bubble nucleation! as well as for 
the droplet nucleation^ is observed in our simulation as 
we are moving along the MFEP of free-energy surface 
F for a single nucleus. Since our MFEP is expected to 
trace the route along the valley on the free-energy surface 
F guided by the phase-field equation Eq. ([1]), we cannot 
answer questions regarding the stability of the MFEP 
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against the thermal fluctuation and the effects of the de- 
pletion of materials to be incorporated into a growing 
nucleus. The latter should be explored using Eqs. ^ or 
dH) instead of Eq. (P). 

The {N, A) cluster model is an attempt to count a fluc- 
tuation of a nucleus around the MFEP that represents a 
more condensed or expanded droplet, for example, in the 
vapor to liquid nucleation. In the (iV, A) two-dimensional 
space, no instability is observed along the A^-axis and the 
instability occurs only along the A-axis when the size A 
is changed while the number of molecule is fixed^. The 
MFEP is the specific route on the (A^, A) two-dimensional 
space on which the size of cluster A is always optimized to 
the given number of molecule N to make the free energy 
minimum. 

In real nucleation phenomena in a conserved system, 
such as the condensation or the cavitation in the liquid- 
vapor nucleation, the evolution of the nucleus does not 
necessarily follow the MFEP on the free-energy surface 
F predicted from the phase-field equation Eq. ^ for the 
non-conserved order parameter. Rather the dynamics 
will be governed by Eq. ^ for the conserved order pa- 
rameter or Eq. ^ of the dynamical DFT when the or- 
der parameter is density. The evolution of the nucleus 
in the conserved system may not follow the MFEP on 
the free-energy surface F but will follow another route 
which also guarantees a monotonically decreasing free en- 
ergy F from the inequality similar to Eq. ([2]). Then the 
criticism that is raised by Lutskoi^ to the {N, A) clus- 
ter model^'^'i for the conserved liquid-vapor nucleation 
system is not correct not only mathematically^ but also 
physically. The {N, A) cluster model may include the es- 
sential feature of real nucleation phenomena governed by 
Eqs. ([3]) or ^ for the conserved system. 

In fact, the real picture of nucleation for the conserved 
system such as the liquid-vapor nucleation will be follow- 
ing: Initially a sub-critical nucleus shrinks or a super- 
critical nucleus grows according to the kinetics governed 
by the non-conserved phase-field equation Eq. ([1]) as the 
depletion of material (monomers) which will be incor- 
porated into the growing nucleus can be neglected. In 
later stage, due to the depletion of material, the dy- 
namics is now switched from the non-conserved equa- 
tion in Eq. ^ to the conserved one in Eq. ^ or to 
the dynamical DFT equation in Eq. For the con- 
served system, the minimum free-energy path (MFEP) 
is merely a virtual path on the free-energy surface F 
that does not take into account the depletion of mate- 
rial due to the increase of the size of nucleus itself or 
the effect of the other growing nuclei within the same 
system. The switching from the non-conserved to the 
conserved dynamics occurs when this depletion becomes 
effective. The distance of this switching point from the 
saddle point of the MFEP depends on many parameters 
of the system considered such as the diffusion length, sys- 
tem size, and the number of nuclei presented etc. The 
study of the MFEP— alone cannot answer this ques- 
tion. This switching point must also be close to the ter- 



mination point of free-energy surface discovered in the 
[N, A) cluster model by Uline and Corti^i^S. Therefore 
the MFEP could be meaningful in the conserved system 
only at the early stage of nucleation for the sub-critical 
nucleus and, probably, only for the super-critical nucleus 
just over the size of the critical nucleus at the saddle 
point. Of course, in the non-conserved system where the 
order parameter represents, for example, the crystallinity 
of the solid-liquid phase transition, the MFEP will be to- 
tally meaningful up to the time when the nuclei start 
to coalescei^. Then the KJMA (Kolmo gorov-Johnson- 
Mehl-Avrami) picture of phase transformationSS will be 
valid in the non-conserved system. 

Finally, it must be noted that the MFEP is usu- 
ful only to characterize the nucleation and the evolu- 
tion of a single nucleus. The nucleation in real ma- 
terials may not simply be characterized by the growth 
of a compact isolated single nucleus. In fact, Shen and 
Debenedetti^ have studied the bubble nucleation in the 
superheated Lennard- Jones fluid, and suggested simulta- 
neous nucleation and coalescence that leads to a ramified 
non-compact nucleus^. Non isothermal effect due to the 
local temperature fluctuation is also sggested to affect the 
nucleation process of bubble recently^. In such a case, 
we have to consider the heat dissipation equation cou- 
pled with the phase field equations, Eqs. ([T|), ([3]) or 
as those equations are considered to describe the isother- 
aml process. 



IV. CONCLUSION 

In conclusion, we have proposed a novel mathemat- 
ical method to explore the minimum free-energy path 
(MFEP) on the free-energy surface of the homogeneous 
nucleation using the phase-field equation for the non- 
conserved order parameter. By employing the local 
square-gradient density functional and by solving the 
phase-field equation for a single circular nucleus using the 
cell-dynamics method, we can calculate the time evolu- 
tion of the size, shape, and free energy of the sub- and 
the super-critical nucleus. This information is used to 
extract the MFEP of the homogeneous nucleation. The 
method is based on the principle in Eq. ^ that the evo- 
lution driven by the phase-field equation for the non- 
conserved order parameter always occurs along the path 
on the free-energy surface along the direction to lower 
the free energy. More rigorous dynamical DFT— or 
the phase-field equation for the conserved order parame- 
ter that can also satisfy an inequality similar to Eq. Q 
may not be effective to prove MFEP though they can be 
used to study the real picture of nucleation and growth. 
Very recently a strategy similar to ours to use the over- 
damped dynamics in the function space was proposed by 
Qiu et al.— where an equation similar to the phase-field 
equation Eq. ([1]) is used as a virtual dynamical equation 
for the discretized free energy to explore the MFEP of 
the nucleation for the capillary condensation. 
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